The possible explanation of electric-field-doped Cqq phenomenology in the framework 

of Eliashberg theory 
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In a recent paper (J.H. Schon, Ch. Kloc, R.C. Haddon and B. Batlogg, Nature 408 (2000) 549) a 
large increase in the superconducting critical temperature was observed in Ceo doped with holes by 
application of a high electric field. We demonstrate that the measured T c versus doping curves can 
be explained by solving the (four) s-wave Eliashberg equations in the case of a finite, non-half-filled 
energy band. In order to reproduce the experimental data, we assume a Coulomb pseudopotential 
depending on the filling in a very simple and plausible way. Reasonable values of the physical 
parameters involved are obtained. The application of the same approach to new experimental data 
(J.H. Schon, Ch. Kloc and B. Batlogg, Science 293 (2001) 2432) on electric field-doped, lattice- 
expanded Ceo single crystals (T c =117 K in the hole-doped case) gives equally good results and sets 
a theoretical limit to the linear increase of T c at the increase of the lattice spacing. 

PACS numbers: 74.70.Wz; 74.20. Fg 
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After obtaining the surface field-effect doping of high- 
quality organic crystals, J.H. Schon, Ch. Kloc, R.C. 
Haddon and B. Batlogg Jjutve carried out the same ex- 
periment on Ceo crystalsclcru, achieving superconductiv- 
ity up to 52 K in hole-doped samplesa. The fundamental 
importance of this experiment is due to the fact that the 
superconductivity is obtained by field-effect doping the 
original single-crystal material and, thus, without modi- 
fying its structural properties. The authors of the exper- 
iments believe as plausible that the doping chasces are 
confined only in a single layer of Cgo moleculcsa. The 
physical parameters of Cgo are rather well known; the 
symmetry of the order parameter appears to be s-wave 
and the electron-phonon interaction is widely indicated 
as the main cause for the appearance of superconduc- 
tivity in this material. In the present paper we show 
that a theory suitable for reproducing the experimental 
T c versus filling 4f^fn|5PPfrt e d in Ref. 1 is the Migdal- 
Eliashberg theoryaB'EHDai3'E3, even though, in fullerenes, 
the Migdal Theorem probably breaks downEJJ due to the 
large energ y o f phonons and the rather low energy of the 
Fermi Ievell3u3. As a first approximation, we neglect 
these complications and see where the theory in its most 
simple form and with the smallest number of free parame- 
ters can lead. Notice that a study of the phenomenology 
of field-doped Cgo in the framework of the Eliashberg 
theory was carried out also in Ref. 14 although in a dif- 
ferent way. Figure 1 (a) was obtained by normalizing the 
original data from Ref. 1 (and from more recent exper- 
imental data in lattice-expanded Cgo, as we will discuss 
later) so that the maximum of T c corresponds to half 
filling (here indicated by n = 1). In order to fit these 
experimental data, we solved the Eliashberg equations 
in the general case of non-half filling, finite bandwidth 
and flat normal density of states. The physical quanti- 
ties that appear in the fit are thus: the bandwidth W, 



the electron-phonon spectral function a 2 F(Q) and the 
Coulomb pseudopotential p. It is also necessary to set a 
cut-off energy u>c related to the renormalization of the 
Coulomb .pseudopotential. We took the value of W from 
literatures : W=250 meV (independent of doping since 
the lattice structure does not changeEl) and we assumed 
the cut-off energy wc=W as in Ref. 10. The spectral 
function is unknown but, as a first approximation, one 
can assume that_aj 2 i 7 '(f2) is proportional to the phonon 
density of statess F(fl) of Ceo (see figure 1 (b)) and 
nearly equal to the transport spectral function a 2 F tr (fl). 

To determine the proportionality factor between 
a 2 F(n) and F(Cl), we took the experimental resistivity 
data p{T) for the electron-doping case from the figure 5 
of Ref. 1, in the range between T = 225 K and T = 245 
K. The p(T) curves appear to have the same shape and 
slope at all gate voltages, thus, as a first approximation, 
the ratio \ tr /ojp between the transport coupling con- 
stant Xtr and the plasma frequency cup can be assumed 
to be doping- independent. The phonon density of states 
F(fl) and the functions a 2 h (Ct) (for both electrons and 
holes) are doping-indepepjdent as well, since the lattice 
parameters do not changeQ. As a result, also the electron- 
phonon coupling constants A e ,h = 2 j dfla 2 h F(fl) /CI are 
independent of doping. From the high-temperature ex- 
pression of the resistivity: p — 87r/csTAt r / (uSps^U 2 ), and 
by using .the value of the plasma frequency reported in 
literatures fiujp — 1.36 eV, we can deduce the valuei-of 
the transport electron-phonon coupling constant At r ll3. 
In the electron-doping case we find Xt r = 1.1 and, as a 
first approximation, we assume Xt r = X e . As far as the 
hole-doped coupling constant Xh is concerned, from resis- 
tivity data and from calculations of the electron density 
of states J.H. Schon et alo estimate that Xh/X e = 1.5 
and so wc find = 1.5A e = 1.65. 

In order to obtain the exact maximum value of the crit- 
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FIG. 1: (a) Experimental T c versus filling n data for hole 
ing (full squares) and electron doping (full circles^fromlil 
(b) The Ceo phonon density of states F(Q) [frornuj]. 



ical temperature in the half-filling case, as shown in fig- 
ure 1 (a), we solved the Eliashberg equations in the half- 
filling case by using for the Coulomb pseudopotential p 
the values 0.395 and 0.336 in the electron and hole doped 
case, respectively. This result, is in very good agreement 
with theoretical calculationsEj that predict p = 0.3 — 0.4. 

Two critical remarks can be made to our basic assump- 
tions. 

First: due to the degeneracy of Ceo-orbitals, one has 
probably to deal with a multi-band problcmEj. It is well 
known that, in such a case, transport and superconduct- 
ing properties can be managed by different groups of 
electrons but, in this way, there would be several free 
parameters and the model would become too complex. 
Instead, our simple approximation allows precise quanti- 
tative predictions. 

Second: we have assumed for simplicity that the nor- 
mal density of states CN.DOS) is flat even if in recent 
theoretical calculationsliZl it is not so. This calculated 
NDOS can be approximately reproduced by a very sim- 
ple analytical formula: Npf(u>) = l+1.5exp(— \u>\/a) with 
a = 10 meV. If we use this expression of the NDOS in 
the Eliashberg equations, with the same values of p, i.e. 
0.395 and 0.336, we find that the exact maximum T c is 
obtained by using A = 0.843 and A = 1.58. The dif- 
ference with respect to the case where Nn(ui) = 1 is 
small, so that the I-V curves measured in the tunable 



SNS weak-link junctions of Ref. 3 can be fitted equally 
well by both the NDOS. Moreover, if the energy depen- 
dence of the NDOS is symmetric and the peak at the 
Fermi level is not too narrow, its effects on some physical 
quantities can be approximately simulated by an efficient 
valueo of A. In conclusion, if AE = (u)peoh 2 )po/7.5 > 
Hlud (where u>n is the Deybe frequency and po is the 
residual resistivity), the effect of a non-flat NDOS can 
be neglectedEl. In our casetl po — 250 ~ 300 pfl cm, 
AE = 62 -~ 74 meV and Tlwd = hu>i n = 27 meV, where 

uj\ og = exp(j J +O ° dfl a In SI). 

Since the lattice is not modified by the field-effect dop- 
ing, we assume the electron-phonon coupling constant A 
and the width of the conduction band W independent 
of filling n. Now we can try to reproduce the exper- 
imental values of T c versus n by solving the Eliashberg 
equations in the non-half-filling caseH. This is a hard task 
because there are three equations to be solved for the cal- 
culation of the gap A(iu> n ), the renormalization function 
Z(iuj n ) and the asymmetric part of the self-energy \{iuj n ) 

- which is always equal to zero in the half-filling case 

- plus one equation that represents the conservation of 
the particles' number and is necessary for calculating the 
shift of the chemical potential 6p. Thus the Eliashberg 
equations in the non-half-filling case, for finite bandwidth 
and s-wave symmetry of the order parameter read:Q: 

+oo 

A(iu) n )Z(iu> n ) = irk B T ^ [X(iu n - iu> m ) - l**/W] ■ 

m— — oc 

fl(\wc\ - u m )P(iu m )Q{iu m ) (1) 

+ 00 

uj n Z(iuj n ) = 7rk B T 2J K iuJ n - ioJ m )N(iuj m )®{iuj m ) 

m=—oo 

(2) 

X(iu n ) = -irk B T 22 H iuj n - iu m )M(iuj m ) (3) 



n = 1 - (irk B T/W) M (^m) (4) 

m— — oo 

where $ is the Heaviside function, iu) m = m{2m + l)k B T 
with m = 0, ±1, ±2, ... and 



M* = /V[l + 



k B Tp 



i?(w m - LU C ) 



(5) 



a^„ - i Wm ) = f°° dn f {n)F{n) (6) 

A(iuj m )Z(iuj m ) uj m Z(iuj m ) 

P(lLO m ) = — ,N(lLJ m ) = — — (7) 



R(iu>„ 



R(iu> m ) = \/u? n Z 2 (iuj m ) + A 2 (iuj rn )Z 2 (iuj rn ) (8) 
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FIG. 2: Coulomb pseudopotential n (open circle), fi* (open 
square) and theoretical fit (solid line) versus filling n: (a) 
electron-doped case, (b) hole-doped case. 

Q(iu) m ) = {arctan(W_(zu; m )) + arctan(W+(io; m ))}/7r 

(9) 



M(iw m ) = 0.5 ln[ 



l + W 2 _{iuj m ) 

i + wi(iLo m y 



(10) 



and W T (iuj m ) = {W T [$V ~ x{iu m )]} / R(iu m ) ■ 

In order to reproduce the experimental values of T c 
versus doping, and having excluded the dependence of A 
on the filling, the only possibility is to assume that the 
Coulomb pseudopotential /i is variable with the filling 
n: fj, = [i(n). We find the values of fJ>(n) by solving 
the equations (l)-(4) conditioned to the obtainment of 
the experimental T c (n). In figure 2 we show the /J.(n) 
dependence in the two cases: (a) electron-doped Cgo and 
(b) hole-doped Cgo (open circles). It is important to 
notice that here the definition of fi* is more general than 
usual because it depends on the shift of the chemical 
potential and on the valuc^of the bandwidth. In the usual 
and less general definitionfl, if lo c = W then fj, = fj,* while 
in our case it is always /i > fi* as we can see in figure 2 
(open squares). 

Is it possible to explain this particular dependence of 
fx on the filling, obtained by the fit of the critical temper- 
ature of figure 1 (a)? This can be .done in two or three 
dimensions. In a very recent paperta the authors affirm 
that it is not so sure that the superconductivity in this 



material is 2D. For completeness we examine both the 2D 
and 3D cases. From the definition of and from the 
very simple analytical expression of the Thomas-Fermi 
dielectric function in the 3D case e(q,n) = 1 + kg(n)/q 2 
it follows: 



H(n) 



1 



2k r . 



e(q,n] 



qdq 



(11) 



where V(q) = 47re 2 /q 2 . In the half-filling case, the cal- 
culation of the integral gives: 



tx{n = 1) 



.k 2 {n 



1) 



8k 2 F 



ln[l + ( 



4fc| 



k 2 (n = 1) 



)] (12) 



where kp and ks are the Fermi and Thomas-Fermi wave 
vector, respectively. We can expand the dielectric func- 
tion in the vicinity of the half filling (n = 1 ) by remem- 
bering that, since the physical parameters A, a 2 F(fl) and 
W are fixed, T c is maximum when /i(n) is minimum, 
which means that e(n) is maximum (see definition of 
/x(n)) i.e. de/dn = for n = 1. For simplicity, and 
to minimize the number of free parameters, we arrest 
the expansion to the second-order terms, even though 
the range of the series expansion is not so small in the 
hole-doping case: 



e{q,n) 



e(q,n= l) + ~ 



1 ,d 2 e(q, n) 



dn 2 



|„=i(n-l) 2 + ... (13) 



and by substituting in the definition of fi(n) we find 



H(ri) = /i(n = 1) 



ln(l + a 2 /[l + b(n-l) 2 }) 
ln(l + a 2 ) 



(14) 



where a = 2kp/kg(n = 1) and b = 
[(l/ks)(d 2 ks/dn 2 )] n =i- By starting from the val- 
ues of fi(n = 1) that give the exact experimental T c in 
the electron- and hole-doping case and using eq. 12 we 
can calculate 2fc^/fcs(« = 1). We obtain 1.05 for hole 
doping and 0.77 for electron doping. Now we are left 
with only one free parameter, b = (1/ ks)[d 2 ks /dn 2 ] n= i, 
that can be adjusted to fit the /j, versus n curves with 
eq. 14. The results are very good (see figure 2 (a) and 
(b), solid lines): for electron doping b = —25 and for 
hole doping b — —4. Of course, here we used the 3D 
version of the Thomas-Fermi theory but, as pointed 
out by the authors of ref. 1-3, the field-effect charge 
injections is likely to be confined to the first atomic layer 
of Cqq crystal, so it is worthwhile to examine also the 
2D case. The analytical expression of the Thomas-Fermi 
dielectric functioned is now e(q,n) — 1 + ks(n)/q and 



k, 



2k- r 



V(q) 



zdq 



(15) 



nhvF Jo e(q, n)y/4k F - q 2 
where V(q) = 2ire 2 /q^. By following the same approach 



as before we find: 

H(n = 1) 



i [arctan( 7 =|= T ) - 1] 



(16) 
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and finally 



H(ri) = /i(n = 1) 



i_ [ arctan ( v#f=r) ~ i] 

t ' [arctai^^l^) - 1] 



(17) 



with a* = a[l + 0.56(n - l) 2 ]. 

In this case, from eq. 16, we can calculate a = 
2kp/ks(n = 1). We find two solutions: 1.9 and 13.8 
for hole doping and 2.1 and 10.3 for electron dop- 
ing. Again we have only one free parameter, b = 
(l/ks)[d 2 ks /dn 2 ] n= i, that can be adjusted to fit the fx 
versus n curves with eq. 17. We find b = —43 and 
b = — 73 for electron doping while for hole doping we 
have b = —4 and b = —8, since now two possible values 
for a are present. In the range of interest, the results 
are exactly equal to the 3D case and the 2D curves per- 
fectly overlap those obtained in that case. If we com- 
pare the 2D and 3D densities of charge carriers obtained 
from the a values of the fits to the values predicted at 
optimal doping and corresponding to 3 carriers per Ceo 
(~ 3 • 10 18 carriers/m 2 and ~ 4.2 • 10 27 carriers/m 3 , re- 
spectively) we obtain interesting results. The 3D density 
and the 2D one obtained from the small a values are of 
the same order of magnitude of the predicted densities, 
even if they don't coincide. On the contrary the 2D den- 
sity obtained from the large a values is more than 2 orders 
of magnitude greater than the predicted one and, thus, 
the corresponding solution for a must be rejected on the 
basis of physical arguments. It follows that the results of 
the fit by using our model do not allow us to understand 
whether the system is better described JpK-SD or 3D ap- 
proach, even if the experimental resultso'Ell clearly point 
to a 2D nature of the charge injected region. 

Note that, surprisingly, the curves T c (n)/T c (n = 1) 
and = l)//i(n) are perfectly superimposed. 

Now, since we have the approximate analytical de- 
pendence of /j from 7i, we can solve the Eliashberg 
equations for any value of the filling and calculate 
different physical quantities. By using the stan- 
dard technique of Pade approximants we can carry 
out the analytical continuation of Z{iuj n ) 1 A(iw rl ), 
xiiwn) and of the normalized superconductive density 
of states (DOS) N s (iuj n )/N N (iuj n ) = [N(iuj n )0(iuj n ) - 
M(iuj n )]/[N]\[(iuj n )<dN(iuj n ) — MN(iu> n )]. Here the quan- 
tities with suffix N are calculated for A(ioj n ) = 0. In 
figure 3, we can see how the shape of Z(u>), A(w) and 
x(^) changes for three different values of the doping in 
the hole-doped case: n = 1 (T c — 51 K), n = 1.189 
(T c = 46 K) and n = 1.287 (T c = 39 K). Figure 4 
shows the normalized superconducting DOS as a function 
of the energy, calculated by the analytical continuation 
of the imaginary-axis DOS for the same doping values 
as in Fig. 3. In the non-half-filling cases (n = 1.189 
and n — 1.287) the calculated DOS is asymmetric. In 
principle low-temperature SIN tunneling experiments on 
this field-doped material would make it possible to ob- 
serve this asymmetry. The calculated values of the ratio 
2A/ksT c are slightly greater than the experimental re- 
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FIG. 3: The real and imaginary part of the functions Z(u>), 
£(w) at T = 6.21 K and of A(lo) at T = 2 K for different 
values of the filling: n = 1 (solid line), n = 1.189 (dash), and 
n = 1.287 (dot) in the hole-doped case. 
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suits of Ref. 3. 

Very recently J.H. Schon et aZ0 have intercalated sin- 
gle crystals of Ceo with CHCI3 and CHBr 3 in order to 
expand the lattice (see figure 1). The maximum of the 
critical temperature (onset of resistive transition) in the 
hole-doped C6o/CHBr3 is 117 K. We can try to explain 




FIG. 5: Coulomb pseudopotential ll (open circle), fi* (open 
square) and theoretical fit (solid line) versus filling n: (a) 
electron-doped case, (b) hole-doped case. 

these new experimental data by following the approach 
previously illustrated. We note that the maximum of 
T c occurs at the same value of doping in the three hole- 
doped and electron-doped cases (see figure 1 (a)). 

As a first approximation we can assume that the 
Coulomb pseudopotential is, for n = 1, the same as in 
the bare Ceo- For simplicity, we also use the electron- 
phonon spectral function of Cgo simply multiplied by a 
constant in order to obtain the correct T c . These two 
approximations can be justified by: (i) the presence of 
an almost constant density of charge carriersaO in the 
bare and intercalated Cgo; (ii) the weak dependence of 
T c from_the details of the actual form of the spectral 
functions (T c strongly depends on LU\ og and on the total 
energy range of the spectral function that are likely not 
to be changed very much by the intercalation) . The val- 
ues of the coupling constant A necessary to get the new 
experimental critical temperatures are 2.11 and 2.76 in 
the hole-doped case (T c = 80 K and T c = 114 K, re- 
spectively) and 1.24 and 1.37 in the electron-doped case 
(T c = 18 K and T c = 26 K, respectively). In fig- 
ure 5 the li and li* values calculated by the solution of 
Eliashberg equations (l)-(4) applied to the T c (n) data 
of lattice-expanded Ceo are shown, as well as the best 
fit obtained by using eq. 14 and eq. 17. Also now the 
fitting curves obtained by using the 2D and 3D formula 
are perfectly superimposed. It is clear that the fits are 
equally good as those shown in Fig. 2. 

Unexpectedly, we find that in all cases the shift of the 
chemical potential at T — T Cl calculated by solving the 
Eliashberg equations, is a linear function of the filling as 
we can see in the six panels of figure 6. The reason of this 
behaviour will be the subject of further investigation. 

Figure 7 (a) shows that the calculated coupling con- 
stant A is a linear function of the lattice constant d. 
We can now extrapolate this behaviour to calculate the 
coupling constant of a hypothetical material with lattice 
constant d w 14. 6A (the increase is of 1% with respect to 
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FIG. 6: The chemical potential shift (open circles), at T = T c , 
calculated from Eliashberg equations and the linear fit (solid 
line), in the hole-doped cases (a,b,c) and electron-doped cases 
(d,e,f). 





Ceo 


Ce,o/CHClz 


Cm/CHBrs 


Ceo/? 


r c m r(K) 


51 


80 


114 


135 




11 


18 


26 


35 


K 


1.65 


2.11 


2.76 


3.51 


Ae 


1.1 


1.24 


1.37 


1.52 


Afc(meV) 


9.9 


16.2 


25.7 


35.5 


A e (meV) 


1.9 


3.2 


4.9 


6.5 


b 3 h D (a= 1.1) 


-4 


-5.5 
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-79 
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d(A) 


14.15 


14.28 


14.43 


14.6 



TABLE I: Important quantities for the six cases examined. 



the case where T c = 117 K)0: \ h = 3.51 and A e = 1.52. 
With these A values and, by using the Eliashberg equa- 
tions, we determine the corresponding critical tempera- 
tures that are 135 K and 35 K, respectively. Finally, fig- 
ure 7 (b) reports the critical temperature as a function 
of the coupling constant, for constant ll and n = 1, in 
the hole- and electron-doped case. The critical temper- 
ature, in the hole-doped case, does not increase linearly 
for A > 3 (d > 14. 5A), in disagreement with the theoret- 
ical predictions of Ref. 12. On the contrary, T c drops to 
zero for A < 0.8 (d < 13.91) as reported in Ref. 20. 

In the inset of figure 7 (b) we can observe that also 
the parameter b is a linear function of the lattice con- 
stant d. All the results are finally summarized in table I. 
In principle, the predictions of Fig. 7(b) and of the last 
column of Table I concerning T c and the gap could be ex- 
perimentally verified in the future. As far as the lattice 
constant dependence of T c shown in Table I is concerned, 
we admit to that the T c enhancement might be not due 
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FIG. 7: (a) (a) Coupling constant A versus lattice constant 
d in the electron-doped case (full squares), hole-doped case 
(full circles) and their linear fits (solid lines); (b) T c versus A 
in the electron-doped case (full squares) and hole-doped case 
(full circles); in the inset, the parameter b is plotted versus d 
in the electron-doped case (full squares) and hole-doped case 
(full circles) together with the linear fits (solid lines). 



to the increase of the lattice constant but due to the en- 
hanced polarization and the resulting enhanced screening 
of the Coulomb interaction. Unfortunately, the formal- 
ism proposed here is not powerful enough to discriminate 
between both scenarios. 



In conclusion, the experimental results concerning the 
complete doping dependence of T c in field-effect hole- 
doped and electron-doped Cgo can be naturally explained 
by the Eliashberg theory generalized to the case where 
the conduction band is finite, non half-filled and the 
Coulomb pseudopotential is filling-dependent in a very 
simple and physically reasonable way. All the values 
of the physical parameters of the present model appear 
plausible. 
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